Deflections of cosmic rays in a random component of the Galactic magnetic field 
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We express the mean square deflections of the ultra-high energy cosmic rays (UHECR) caused by 
the random component of the Galactic magnetic field (GMF) in terms of the GMF power spectrum. 
We use recent measurements of the GMF spectra in several sky patches to estimate the deflections 
quantitatively. We find that deflections due to the random field constitute 0.03 — 0.3 of the deflections 
which are due to the regular component and depend on the direction on the sky. They are small 
enough not to preclude the identification of UHECR sources, but large enough to be detected in the 
new generation of UHECR experiments. 



I. INTRODUCTION 

The Galactic magnetic field (GMF) plays an important 
role in the propagation of cosmic rays even at highest en- 
ergies. Expected deflections — of order few degrees or 
larger — are comparable or exceed the angular resolu- 
tion of the existing cosmic ray experiments. Such deflec- 
tions may therefore be observable. Their understanding 
is crucial when searching for sources of the highest-energy 
cosmic rays if the latter are charged particles. 

The detailed study of deflections of ultra-high energy 
proton primaries in the GMF is particularly important if 
the deflections in extra- galactic magnetic fields are small. 
According to the results of Refs. p| this is likely to be 
the case (see, however, Ref. @). 

The Galactic magnetic field has been shown to have 
both regular and turbulent components. The regular 
component is thought to have a spiral structure remi- 
niscent of the Galactic arms with one or more reversals 
toward inner (and probably also outer) Galaxy and the 
magnitude of order 3 fiG in the vicinity of the Earth Q . 
The corresponding global model of GMF was constructed 
in Ref. |^ and is often used in the discussion of propa- 
gation of cosmic ray primaries in the Galaxy. According 
to this model, protons with energy 4 x 10 19 eV can be 
deflected in the regular GMF by ~ 5°. There are indi- 
cations that such a coherent deflections may indeed be 
present [E IE HI in the cosmic ray data. 

The random component of GMF causes the spread of 
arrival directions of UHECR around the mean position, 
thus diluting (and potentially destroying) important in- 
formation about the actual location of the source [jj . Un- 
der certain conditions on the magnetic field it may also 
lead to the "lensing" of cosmic rays [ToL HTI| provided the 
number of sources contributing to the observed UHECR 
flux is small, as is favored by the statistics of clustering 
|12| . Thus, the influence of a random fields is not sim- 
ply destructive, but may give useful information on the 
cosmic rays and GMF itself. 

Observationally, the magnitude of the random com- 
ponent of GMF is comparable to the magnitude of the 
regular one. However, the deflections of cosmic ray pri- 



maries in the random field are expected to be consider- 
ably smaller. Indeed, if the correlation length L c of the 
random component is much smaller than the propaga- 
tion distance D, the deflections caused by the random 
field are proportional to \JDL Cl see e.g. Ref. [T^|. while 
the deflections in the regular field are proportional to the 
distance D itself. Usually the deflections in the random 
field are estimated as (see e.g. Ref. |l4j ) 
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where E is the energy of CR primary and B r is the rms 
value of the random magnetic field strength. However, 
the global picture of GMF is very complicated. The mag- 
netic field may be structured on relatively large scales, in- 
cluding possibility of magnetic bubbles, sheets, filaments, 
magnetic "winds" etc., so one cannot rely on (and be re- 
stricted to) a generic estimate uniformly over the whole 
sky. 

The rms deflections of cosmic ray primaries in the ran- 
dom field are determined unambiguously by the magnetic 
field power spectrum. The latter can be extracted from 
observations, see e.g. Ref. T5j. Several existing mea- 
surements of the GMF power spectrum indicate that the 
correlation length can be large in some directions on the 
sky. In this case the deflections are no longer propor- 
tional to \JDL C , and the estimate Q does not apply. 

The purpose of the present paper is to estimate the 
deflections of UHECR primaries in the turbulent compo- 
nent of the Galactic magnetic field taking into account 
possible dependence of GMF parameters on the region of 
the sky and the possibility of large corelation length. We 
express the UHECR deflections directly through the ob- 
servable parameters characterizing the power spectrum 
of the magnetic field. To this end we derive the relation 
between the mean square deflection and the power spec- 
trum (PWS) of the GMF fluctuations. The result is most 
conveniently represented through the factor R defined as 
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where 5 r and 8 U are deflections in the random and uni- 
form components of GMF, respectively, and B u ± is a 
projection of a uniform field onto the direction orthogo- 
nal to the line of sight. The factor R varies between 
and 1 and is expressed in terms of the power spectrum of 
the random field by Eq. I)28[l . Moreover, we show that it 
is a function of a single observable parameter 8 Cl the an- 
gular scale of the break in the relevant structure function 
of GMF. In the case of small correlation length when the 
estimate applies, one has R ~ \fW c . The values of 
R in our Galaxy derived from the existing observational 
data vary in the range R ~ 0.03 — 0.3. This implies typi- 
cal deflections of a 4 x 10 19 eV proton in the random field 
of order 0.2° — 1.5° depending on the direction. 

The paper is organized as follows. In Sect. Ill Al we 
recall the relations between the power spectrum and the 
correlation length and introduce the notations. In Sec. lIIII 
we describe existing observations of the random magnetic 
field, in particular, the random-to-uniform ratio and the 
parameters of the power spectrum. In Sect. II VI we turn to 
the deflections of UHECR in the random magnetic field 
and derive the expression for the coefficient R. Sect. Ivl 
summarizes our results. 



II. TURBULENT FIELD 



correlation function which can be written as 



(B a (q)B* b (q')) = ^-(S ab 



n a n b )8 3 (q-q'), (3) 



where n a = q a /q is a unit vector in the direction of q and 
the projection tensor ensures the divergence-free nature 
of the magnetic field, q a B a (q) = 0. The dimensional 
normalization factors are chosen in such a way that the 
power spectrum B(q) has physical units of B 2 , i.e. in 
our case it is measured in the units of (Gauss) 2 . The 
correlation function of the magnetic field fluctuations is 
defined as 



£(r) = (B a {v )B a {v + r)> = J ^ B(q) e"*' 
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It determines the rms value of the field amplitude B T1 



B 2 r = (B a B a ) = £(0) , 



and the correlation length L c , 
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A. Correlation length and power spectrum of GMF 

To introduce notations and specify our assumptions, 
let us first consider the statistical properties of a random 
magnetic field B a (r), where a = 1,2,3. We define the 
Fourier components of the magnetic field according to 



B.(r) = J d 3 qB a (<i)< 
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Here B a refers to the fluctuating component of the to- 
tal magnetic field; the regular part of GMF has to be 
treated separately. In what follows we assume that the 
fluctuations obey Gaussian statistics and are spatially 
homogeneous and isotropic. 

The last two assumptions deserve a comment. The 
statistical characteristics of the magnetic field are differ- 
ent in different sky patches. Our analysis and results 
should be applied to each of these patches separately. 
We assume that statistical properties of GMF are (ap- 
proximately) constant over a single patch. Also, within 
one patch the magnetic field fluctuations may not be 
isotropic, with the preferred direction being set by the 
regular component of GMF. The present data are not 
sufficient to establish or rule out the isotropy of GMF 
fluctuations. With the more precise data this assump- 
tion may need to be reconsidered, and the analysis may 
need to be refined. 

With the above assumptions, all correlators of the 
magnetic field can be expressed in terms of the two-point 



The energy density contained in a random component 
of the magnetic field is related to the field variance as 
p B = B 2 /8n. 

As will be discussed in Sect. lIII 51 in a certain range of 
momenta the existing observations support a power-law 
behavior of power spectra of the magnetic field fluctua- 
tions, 



1 



(7) 



Since the energy density in the magnetic field is finite, 
there has to be a break in the pure power-law behavior, 
which can be parametrized as 



B(q) = 



Al q -1 



OL\ — 1 



at q < q c 



at q > q c 
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where A is a normalization constant and q c is the mo- 
mentum scale at which the break in the spectrum occurs. 
(Note that abrupt ultraviolet and infrared cut-off can be 
modeled as ct2 — ► oo and a% — > — oo, respectively.) In 
what follows we assume this form of GMF power spec- 
trum. The variance B r , and consequently the energy 
density, converges if a>i < 1 and a.^ > 1. Summing up 
contributions from both parts of the spectrum one finds 
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FIG. 1: Schematic view of the turbulent field and the under- 
lying geometry. At small angles, 6 < C = Le/D, there are 
correlated fluctuations along the entire length of the lines of 
sight. At large angles only neighboring structures contribute 
to correlations. 

Finiteness of the correlation length requires stronger con- 
straint, ot\ < 0. One then has 

Lc=JL {ai -l) { a 2 -l) 
2q c a\a 2 

If < ot\ < 1, the correlation length diverges at small 
momenta and is dominated by the largest possible dis- 
tance scale in the problem. Deflections of cosmic rays 
are most significant in this case. 

B. A physical picture and observables 

A toy model leading to the PWS with the break is 
illustrated in Fig. Here the turbulent magnetic field 
is contained within the bubble-like structures. Denote 
by Le the typical size of the largest energy-containing 
eddies, which is also called an outer scale of turbulence 
or energy injection scale. It is often assumed that Le 
corresponds to a typical size of the supernova remnants, 
and that inside eddies the PWS corresponds to MHD 
or Kolmogorov turbulence. With the assumption that 
fluctuations in different "bubbles" are uncorrelated, the 
corresponding correlation function can be approximated 
as 

£(r) = R 2 J 3 - r 2 / 3 at r < R E , (u) 
f(r) =0 at r > R E , y ' 

where Re = Le/2. Using Eq. JJJJ we see that the corre- 
lation length of the magnetic field fluctuations is related 
to the scale Le as L c = Le/5. Similar expression for the 
correlation function can be also obtained from Eq. 
assuming Kolmogorov spectrum, oii = 5/3, with infrared 
cut-off at a momentum scale q c ~ 2tt/Le- 

In this example the 3d correlation function is zero at 
large separations. In general, however, the fluctuations 
may be present on large scales as well. This is the case, 
e.g., when the Kolmogorov turbulence is contained not 
within the "bubbles" as above, but within long sheets or 



filaments. When substantial power is present in fluctua- 
tions up to the largest scales, the correlation length may 
diverge. 

Three-dimensional power spectrum of the magnetic 
field fluctuations, B(q), and the corresponding 3d cor- 
relation function are not measured directly. Instead, one 
measures the parameters of the two-dimensional angu- 
lar correlation function K(9) of some physical observ- 
able which can be, e.g., the intensity of the polarized 
synchrotron radiation or Faraday rotation measure along 
different directions, see Sect. Hill for examples. For such 
observables the relation between 3d correlation function 
and the angular correlation function involves integration 
along the line of sight as illustrated in Fig. ^ For ex- 
ample, in the case of Faraday rotation measure of extra- 
galactic sources such relation reads 

K(9)oc /^dndra^Mdn-ral). (12) 
Jo 

The directly measurable quantities are the structure 
function 

±S(6) = K(0)-K(9), (13) 

and the two-dimensional power spectrum Ci which at 
small angles is related to K(9) by the Hankel transform 

K{e) = ±- / d 3 {19) Idl. (14) 
27T Jo 

Here the multipole I corresponds to a typical angular 
scale of 9 — ir/l . 

In the case of the power-law behavior of the power 
spectrum, Eq. J7J), both S(9) and Cj also follow power- 
laws, S(9) cx 9° and Ci cx I 1 . In a certain range of 
a (see Appendix El for details) the three exponents are 
related as 

13 = a, (15) 
7 = a + 2. (16) 

Eq. I|15|l is valid when < a < 2. When the correla- 
tion length is finite, K(9) — > at sufficiently large an- 
gles, and the structure function approaches constant (cf. 
Eq. (|13|) h Therefore, (3\ = for any negative a±. On 
the other hand, at negative a the relation i|l(j|) holds. 
Making use of these two relations, the exponent a can be 
deduced from the observable quantities f3 and 7. How- 
ever, systematic observational effects (e.g., Faraday depo- 
larization and finite beam width in the case of polarized 
synchrotron radiation) should be absent or accounted for. 

The momentum scale q c at which the break in B(q) 
occurs is not observed directly either. Instead, one ob- 
serves the break in the power-law behavior of S(9) at 
some angular scale 9 C (or the break of C; at some multi- 
pole l c ). The scale q c can be estimated using the relation 
l c = Dq c /2, or 

2tt , s 

= oj> ' (17) 
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see Appendix [X] These relations involve another un- 
known parameter, the propagation distance in the mag- 
netic field, D. Fortunately, in the expression for the cos- 
mic ray deflections q c and D enter as the product q c D, 
so the result can be conveniently expressed in terms of 
the observable parameter 9 C . 

In the model of Fig. ^ the structure function at small 
angular scales 9 < 6 C = Le/D should reflect Kolmogorov 
turbulence, 02 = 5/3. On larger angular scales the struc- 
ture function should gradually become flat, (3\ — > 0. The 
correlation function in this case decays on large scales 
as K(6) cx which corresponds to Ci oc Z _1 |16|. In 
other words, existence of small (compared to D) correla- 
tion length corresponds to fli = and 71 = 1. 



III. OBSERVATIONS OF RANDOM 
COMPONENT OF GMF 

Current knowledge of the Galactic magnetic field is 
based on: (i) Faraday rotation measurements of Galac- 
tic and extragalactic radio sources, (ii) starlight polar- 
ization data, and (iii) observations of diffuse Galactic 
synchrotron emission. Different methods are sensitive 
to the magnetic field in regions with different physical 
conditions. Faraday rotation is sensitive to a field in a 
warm ionized medium, stellar polarization measurements 
sample the field in regions occupied by interstellar dust 
grains, while synchrotron radiation originates from re- 
gions containing fast electrons. 

Faraday rotation measure (RM) is sensitive to the pro- 
jection of the magnetic field on the line of sight. The field 
direction along the line of sight is given by the sign of RM. 
The magnitude of the random field B r can be estimated 
by analyzing the deviations of RM from the uniform field 
along different directions. 

Stellar and synchrotron polarization data contain in- 
formation about the field perpendicular to the line of 
sight. This property is convenient for our purposes since 
the plane-of-the-sky component of the magnetic field de- 
termines also the deflections of UHECR primary par- 
ticles. The ratio of the amplitudes of the random to 
the uniform magnetic field components can be estimated 
along a single direction. To extract the power spectrum 
one has to study angular correlation function of the po- 
larization data. 



A. The relative strength of uniform and random 
fields 

Synchrotron emission. The total magnetic field 
strength is related to the synchrotron emissivity, while 
the polarization of the Galactic diffuse synchrotron back- 
ground offers a method for determining the ratio of uni- 
form to random field strengths. Namely, the observed 
fractional polarization p Q b s along a given direction obeys 



m 

£°* = ^ , (18) 

Pmax B uX + B r± 

where the subscript _L on B u and B r refers to the plane- 
of-the-sky components of uniform and random field, re- 
spectively. In Eq. (|18|l p max is the fractional polariza- 
tion that would be observed for a perfectly uniform field, 
Pmax ~ 0.72. The fractional polarization of p Q b s ~ 35% 
was found in Ref. |18| to be a typical maximum for our 
Galaxy. This implies 

B u ,±/B r , ± « 1 . (19) 

Note that this ratio depends upon direction. For in- 
stance, in the direction of the Galactic anti-center p b s ~ 
20% (being averaged over -20° < b < 20°), which cor- 
responds to B u ,x/B r ,±. ~ 0.62. The typical coherence 
length was estimated in Ref. [18| to be less than 75 pc, 
while the distance to the region where polarized emission 
originates was found to be about ~ 0.5 kpc. 

Starlight polarization. Polarization in starlight appears 
because of selective absorption by interstellar dust grains 
whose minor axis is aligned with the magnetic field B. 
The same expression, Eq. (|18|l . is valid for the starlight 
polarization data as well (with p max being related to dust 
extinction) . The resulting magnitude of the random com- 
ponent of magnetic field derived from the starlight polar- 
ization data is consistent with Eq. i|19|) . For example, the 
estimate of Ref. |l]j reads B U> ±J B rt ± ps 0.8. 

Faraday rotation. Unlike polarization data, the Fara- 
day rotation measure is sensitive to the magnetic field 
component parallel to the line of sight. Another disad- 
vantage of this method is that it does not allow to find 
the ratio of random to uniform components of the mag- 
netic field along a given direction. However, this infor- 
mation can be extracted from the residuals of a fit to a 
uniform field provided RMs in many neighboring direc- 
tions are known. For instance, such kind of study for a 
particular region of the sky of about 10° x 10° centered 
at (l,b) ps (140°, -40°) was carried out in Ref. [3. It 
was found that B r ss B u , in agreement with Eq. iJTSj l. 

B. The power spectrum of magnetic field 
fluctuations 

Synchrotron emission. The statistical properties of po- 
larized synchrotron emission depend upon direction on 
the sky and are different for different observables. In 
Ref. [2(| the angular power spectra (APS) of the Parkes 
survey of the Southern Galactic plane at 2.4 GHz were 
analyzed. It was found that in the multipole range 
40 < I < 250 (0.7° < 6 < 5°) the APS of E and B com- 
ponents of the polarized signal has the slope 7 ps 1.5, and 
the power spectrum of polarization angle corresponds to 
7 ps 1.7. Similar results were found in Refs. |2l], [23 for 
other Galactic latitudes. In particular, while being close 
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to 1.5 on average, the slopes of E and B components in 
the multipole range I < 1000 (9 > 10') were found |2l| to 
be in the range 1 < 7 < 2.7 depending on the particular 
region of the sky and the survey used. 

Negative values of a (derived with the use of Ea. (jlfifl . 
if applicable) indicate that in many sky patches the cor- 
relation length of magnetic field may be small. 

Starlight polarization. The angular power spectrum 
of the starlight polarization for the Galactic plane data 
< 10°) is consistent with 7 « 1.5 for all angular scales 
9 > 10' (or I < 1000), see Ref. [H. 

Faraday rotation. Structure functions of the rotation 
measure of extragalactic radio sources were studied in 
Refs. [l5l 123 . l24l l25l 1291 . Three different sky patches were 
considered in Ref. [23 ■ In two patches the index j3\ of the 
structure function of rotation measure was found to be 
consistent with zero on large angular scales > 2°, while 
in the third positive (3 was observed. There is a clear 
drop in the structure function on small angular scales 
9 < 0.1°, [HE]]]. Therefore, the break in the spectrum 
has to be at angular scales 0.1° < 9 C < 2° (9 C cannot be 
quantified more precisely as there are no data points at 
these intermediate angular scales). 

In Ref. [2^| shallow (with fli < 0.3) structure func- 
tions in several sky regions near the Galactic plain were 
found over the range 0.3° < 9 C < 10°. Similar result was 
obtained in Ref. J25j. 

Minter and Spangler have studied structure func- 
tions of RM from polarized extragalactic sources for a 
particular region with previously mapped emission mea- 
sure of warm ionized medium. In this paper fluctuations 
in electron density were factored out and the power spec- 
trum of fluctuating magnetic field was determined. The 
spectrum of random magnetic field derived in Ref. fl5| 
can be parametrized by Eq. (jHJ with A ss 4.5 x 10 -2 \xG 2 . 
At large scales the angular structure function is con- 
sistent with the two-dimensional turbulence, ol\ = 2/3, 
while at small scales q > q c the spectrum coincides with 
the Kolmogorov turbulence U2 = 5/3. The break in 
the spectrum occurs at 9 C ~ 0.07° which corresponds 
to 2ir/q c 3.6 pc assuming D = 3 kpc. With the pa- 
rameters found in Ref. [lj|, Eq- © gives B r « 1.6 fiG. 
Note that for the uniform component of the magnetic 
field in the same region one has B u 2.2 /iG and 
B u .\\ ~ —0.8 \iG. The slope of a\ = 2/3 was measured 
up to 2ir/q ~ 80 pc. Thus, in this particular sky patch 
the correlation length of magnetic field fluctuations either 
diverges or is larger than 80 pc. 

Power spectra and the structure functions of the ro- 
tation measure of the diffuse Galactic polarized radio 
background in several sky patches near the Galactic plain 
were studied in Ref. |27| . Angular power spectra show a 
spectral index 71 ps 1, while the structure functions are 
approximately flat, /3i « in the range 0.1° < 9 < 10°. 
This is indicative of the field uncorrelated on large scales. 
The structure functions may show a break at 9 C close to 
0.3°, which is at the same spatial scales (« 3.9 pc assum- 
ing a path length of 600 pc) as a break in the structure 



function in the RMs of extragalactic sources of Ref. [lfj . 
Note that the RM of extragalactic sources probes the 
complete line of sight through the Galaxy, whereas, as 
a result of depolarization, the synchrotron emission ob- 
served at low frequencies only probes the nearby interstel- 
lar medium. In a similar study |28l | in anther sky region, 
a shallow structure function of the rotation measures was 
found, (3\ ~ 0.2 for 4' < 9 < 5°, while at smaller angular 
scales the structure function steepens. 

IV. UHECR DEFLECTIONS IN THE RANDOM 
MAGNETIC FIELD 

In this section we show that the knowledge of the ra- 
tio B u ^±/B r , the exponents ol\ and a2, and the angular 
scale 9 C is sufficient to quantify the spread of deflections 
of UHECR primaries caused by the random component of 
GMF. As we have seen in Sect. IIII Al the existing obser- 
vations suggest that the magnitudes of the random and 
uniform components of GMF are comparable. In what 
follows it will be convenient to normalize the deflection 
due to the random field to the deflection S u which would 
occur in the uniform field over the same distance and at 
the same particle energy and charge. After traveling the 
distance D in a uniform magnetic field, a particle with 
the electric charge Ze and energy E is deflected by an 
angle [3j| 

5 U = ^ B u>x . (20) 

This has to be compared to the mean square deflection 
angle 8 r in the random component of the Galactic mag- 
netic filed. 

A. Mean square deflection in terms of magnetic 
field power spectrum 

Deflections of UHECR primaries by random magnetic 
field were studied in many papers, see e.g. Refs. |M 
I30L l3ll I32I l33l 133 . 135L l36j . However, usually the main 
focus is the diffusive regime in the extra-galactic mag- 
netic field. Deflections in the turbulent component of 
GMF were studied in Refs. [ToL Hl|. with the emphasis 
on the possibility of magnetic field reconstruction with 
future high statistics cosmic ray data. A generic turbu- 
lent component of GMF with a simplifying assumption of 
a cell-like structure was also included in the Monte-Carlo 
simulations of Ref. 37]. To our knowledge, estimates of 
UHECR deflections in the random field based on mea- 
surements of the MF power spectra in specific sky patches 
do not exist in the literature. The UHECR deflections 
in the situation when the coherence length is not small 
(which might be relevant for the case of realistic GMF) 
was not studied in detail either. 

Propagation of UHECR primaries in the Galaxy is 
quasi-rectilinear, with typical deflection angles not ex- 
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ceeding 10° — 20° even for lowest energies. The contri- 
bution of turbulent field in these deflections is expected 
to be even smaller. Therefore, a ballistic approximation 
gives a good description of UHECR propagation. In this 
regime, the deflection angles are characterized by the fol- 
lowing line integrals, 



5 ^ C 
i ~ ~E 



D 



dz e ik B k (z) , 



(21) 



where the axis z is chosen along the particle trajectory 
and indices i, k — 1,2 label two orthogonal directions. 
The mean square deflections are 



Z 2 e 2 



dz dz' e ik e jp (B k (z)B p (z')) 



(22) 

Here the average is taken over the ensemble of different 
realizations of the turbulent magnetic field B a (x) . For a 
statistically homogeneous random field the correlator in 
Eq. 122|) is the function of r = z' — z 

e ik e jp {B k {z)B p {z')) = £ u (r) + £ 22 (r) = £±(r) , (23) 

where £a(r) = (Bi(z)Bi(z + r)) (no summation over i). 
Using Eq. Q which enforces the divergence-free con- 
straint one finds 



£i(r) = 4tt 



dq 



10 Q 
This relation implies 



B(?) 



sin(qr) cos(gr) sin(gr) 



qr 



1 



q 2 r 2 



q3 r 3 



dr £j_(r) = - dr f(r) 



(24) 



Note that the assumption of chaotically oriented mag- 
netic cells, which is often made, would give instead 
£±( r ) — (2/3) £(r). However, this assumption is incon- 
sistent with the divergence-free nature of the magnetic 
field. 

:' to r and 



Changing variables in Eq. I|22ll from 
= (z + z') one obtains 



2Z 2 p 2 
2 l/j e 2 

d r - B r 



„ du L dr W> (25) 



where Eq. |J5} was used. It is convenient to represent 
the result as a ratio of rms deflections in random field to 
the deflection in the uniform field S u given by Eq. (|20|l . 
Thus, we arrive at Eq. J2J) where the dimensionless factor 
R is 



R 2 



2 

D 2 



„ du h dr «0) 



(26) 



This factor varies between zero and one. 

If the correlation length L c defined by Eq. JJjJ is much 
smaller than the propagation distance D, the upper limit 



in the integral over r can be extended to infinity. One 
then finds 



R 



D 



(27) 



In the general case, the expression Eq. H26J1 can be 
brought to the form 



R 2 



4tt 



dq 



B(q) f(Dq) , (28) 



where 



f{x) = Si(x) + 



cos x sin x 



(29) 



and Si(x) = L dy sm(y)/y is the integral sine function. 
At small arguments the function f(x) grows linearly as 
f(x) = 2x/3+0(x 3 ), while at x }Z 2tt it rapidly converges 
to the asymptotic value n/2. 



B. Mean square deflections in the random 
component of GMF 

With the assumption that the power-law spectrum of 
the turbulent magnetic field is given by Eq. (JSJ, as sup- 
ported by the existing observations, Eq. I|28() gives 



r2 _ (gg - 1)(1 - ai) x 



(a 2 - ax) 



(30) 



(Dq c ) ai 



CK1 — 1 



Dq c 



dyf(y) 

,A+ai 



(Dq c ) c 



dyfiy) 



Dq c y 



l+a 2 



where f(y) is defined in Eq. (|29|) . As one can see, the 
final result depends on the product Dq c . Therefore, with 
the use of Eq. I|17fl it can be rewritten in terms of the 
single directly observable scale 9 C — 2ir/Dq c . 

In the case Dq c ^> ir and a\ < (when both integrals 
are saturated at y = Dq c ), we recover Eq. (|27fl. 



R 




9 C (qi - l)(q 2 - 1) 
4 ai(X2 



(31) 



For ot\ varying within 0.2 < ot\ < 0.8 and a 2 = 5/3 
(which corresponds to the Kolmogorov turbulence) , there 
exists an approximate analytic expression for R which 
holds with an accuracy of about 10%, 



R « (Dq c ) 



(«i-l)/2 



(9c/27T 



,(l-ai)/2 



(32) 



In the general case the factor R has to be calculated nu- 
merically. Its dependence on a,\ for a 2 = 5/3 in three 
cases 6 C = 6°, 8 C = 0.6° and 8 C = 0.06° is shown in Fig.H 
by the dotted, dashed and solid lines, respectively. We 
recall now that in many "low" resolution all-sky studies 
the steepening of APS is not detected, up to large mul- 
tipoles, I ~ 1000. This suggests that C < 10' and the 
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1 



R 




Ooi l 1 1 1 1 1 1 1 1 l_ 

-0.8 -0.4 0.4 0.8 

OCj 



FIG. 2: The coefficient R as a function of ai is shown by 
dotted, dashed and solid curves for 8 C = 6°, 0.6° and 0.06° 
respectively. The data-points correspond to the APS derived 
for sky regions A-E, as discussed in the text. These regions 
are displayed in Fig. |3] 



dashed line in Fig. [21 may serve as an upper limit for the 
factor R. 

The labeled data-points in Fig. |3 represent the data 
based on the Faraday rotation measurements. The cor- 
responding survey regions are shown in Fig.|3|as colored 
patches labeled in the same wa y a s the data-points by 
A (Ref. Ill), B1,B2 (Refs. |H, C1,C2 (Ref. H3), 
D (Ref. |H) and E1,E2,E2 (Refill). Only in the re- 
gion A the transition to the Kolgomogorov turbulence 
was detected at the scale 8 C ~ 0.07°. Though in some 
other regions an indication for the break was observed, 
the transition to the spectra with ot2 > 1 was not es- 
tablished. The corresponding limits on 9 C are shown as 
downward arrows. As explained in the Appendix ^ any 
negative a corresponds to j3 — in the structure func- 
tions. Therefore, we plot data-points with (3 = as down- 
ward arrows at ct\ = turning to the left (pointing to the 
region of negative a\). Finally, in the regions C1,C2 both 
the structure functions and the APS of rotation measure 
were obtained |23 resulting in (3 w and 7 w 1. This 
allows to specify a\ in these regions as a\ & — 1. 

Small observed values of (or upper limits on) 9 C indi- 
cate that either the scale 2w/q c is small, or the extent of 
GMF along given direction, D, is large. In either case 
the resulting coefficient R is small, 0.02 < R < 0.2. Note 
that the application of Eq. to the directions along 
the Galactic plane should be done with care. Namely, 
deflections in the regular field cannot be approximated 
by a simple relation J5DJ), but should be replaced by the 
integral along the line of sight. Existing field reversals 
can diminish deflections in the regular field considerably. 



V. CONCLUSIONS 

Deflections of UHECR in the random component of 
the Galactic magnetic field are usually discussed in the 




FIG. 3: Regions A-E in Galactic coordinates where the break 
in the APS was detected. Corresponding values of R for these 
regions are shown as data-points in Fig. [5] 



limit when the correlation length L c is much smaller than 
the propagation distance D. However, even in this limit 
both the correlation length and the propagation distance 
in GMF (and therefore typical deflections) vary with di- 
rection. Moreover, the existing GMF data suggest that 
the assumption of small correlation length may not be 
valid uniformly all over the sky. 

We have calculated deflections in a more general ap- 
proach which does not require correlation length to be 
small and relies directly on the spectrum of GMF fluc- 
tuations measured in a relatively small patch of the sky. 
Our method therefore automatically takes into account 
variations of the GMF characteristics with direction. We 
have shown that the ratio of the deflections in the random 
and uniform components of GMF, Eq. @, is expressed in 
terms of the factor R which depends on the spectrum of 
the magnetic field fluctuations as given by Eq. I|28|l . For 
the power-law spectrum with a single break the factor R 
can be written as a function of one directly observable 
parameter 9 C , the angular scale of the break in the rele- 
vant structure function. In the case of a small correlation 
length one finds R ~ \[W C . 

Using the measurements of the GMF power spectrum 
in the sky regions where it is available, we have shown 
that the deflections in the random component are small, 
0.03—0.3 of the deflections in the uniform field, see Fig- El 
This is sufficiently small not to preclude identification 
of sources of UHECR using methods of Refs. 0, S H|. 
For instance, the deflection of a proton with energy 
E = 4 x 10 19 eV due to the random component of GMF is 
expected to be about 0.2° — 1.5° depending on the direc- 
tion. This is below the resolution of the AGASA experi- 
ment, but can be above the resolution of the HiRes detec- 
tor in the stereo mode and the expected resolution of the 
Pierre Auger experiment. Thus, the detailed study of the 
random component of GMF is particularly important for 
the interpretation of data which will be collected by the 
new generation of UHECR experiments. To this end, the 
all-sky map of the essential parameters determining the 
power spectra of GMF is highly desirable. These maps 
may be obtained from the measurements of Faraday rota- 
tion and maps of diffuse polarized synchrotron Galactic 
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emission. 
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APPENDIX A: ANGULAR STRUCTURE 
FUNCTIONS, ANGULAR POWER SPECTRUM 
AND UNDERLYING 3D POWER SPECTRUM. 

In this section we derive Eqs. ((15(1 - 1(17(1 . These equa- 
tions are valid for a large class of observables including 
Faraday rotation measures. To simplify the presentation 
we derive them in the case when the correlation function 
is given by equation Q. 

The relation between 3d correlation function and an- 
gular correlation function involves the integration along 
the line of sight, 



K(0) 



D 



r1dr ir ^dr 2 £(|r x - r 2 | 



(Al) 



where 9 is the angle between the vectors ri and r 2 and a 
takes the values and 2 in the case of point sources and 
diffuse radiation, respectively. We concentrate in what 
follows on the case a = relevant for all points in Fig. |3 
except Ci t 2 and D. 

Introducing the variables v = r± — r 2 and y = (ri + 
r 2 )/2 and taking the limit of small angles one can write 
Eq. HAlfl as follows, 

fD/2 r 2y 

K{6) = 4 dy dv £(vV+|/ 2 6> 2 ). (A2) 

Jo Jo 



Integration over dv gives 



fD/2 , 

K(0) = 16tt / dy 

Jo Jo 



7T 7T 

2 Mvtf) - 



Si(2yq) 



(A3) 



where Jo and Si are the Bessel and integral sine func- 
tions, respectively. According to definition 1)13(1 . this re- 
lation gives for the structure function S(9): 



fD/2 , 

S(0) = 16tt 2 / dy 

Jo Jo 



d 4B( q) [1 



MyqO)]- (A4) 



Let us assume that B is a broken power law, Eq. J5J|. 
Then the structure function is also a broken power law 
behaving as S cx Q^ 1 and S oc 0^ 2 below and above the 
angular scale 8 C , respectively. The exponents 0i are re- 
lated to exponents ai in Eq. JHJ as follows, 



= 2 if a > 2, 

= a if < a < 2, 

0=0 if a < 0. 



The break occurs at 



2tt 



(A5) 



(A6) 



The angular power spectrum, Ci , and the angular cor- 
relation function at small angles are related by the Hankcl 
transform, Eq. 114(1. If B is given by the power law © 
with ol\ < (i.e., the correlation length converges), the 
integral in Eq. (|A3|) is saturated in the region yq 3> 1 . In 
this region Si(2yq) — > ir/2; therefore, the last two terms 
in Eq. l|A3j) can be neglected. (This holds whenever the 
correlation length is much smaller D and the integration 
over dv in Eq. I(A2(I can be extended to infinity.) Com- 
parison of Eqs. I(A3(I and l(14|) gives then 



C t = 



16tt j 



D/2 n 

b(L 

o \y 



ydy. 



(A7) 



Note that Eq. I(A7(I is generalized trivially to the case of 
non-zero a by replacing ydy by y 1+2a dy. For a broken 
power-law behavior of B, the angular power spectrum is 
also a broken power law, Ci oc l^ 1 , where the exponents 
7i are 



7i = a i + 2 
and the break occurs at l c = Dq c /2. 



(A8) 
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